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Possibility of the structured mixed phases at first order phase transitions in neu- 
tron stars is reexamined by taking into account the charge screening effect. The 
Maxwell construction is shown to be not conflicted with the Gibbs conditions once 
the Coulomb potential is properly taken into account. Taking the hadron-quark 
deconfinement transition as an example, we explicitly demonstrate a mechanical 
instability of the geometrical structure of the structured mixed phase by the charge 
screening effect. In this case we have effectively the picture given by the Maxwell 
construction. 



X 



1. Introduction 

It is now commonly accepted that various phase transitions may occur in 
compact star interiors or during the gravitational collapse from progeni- 
tor stars. Possibilities of the meson (pion and kaon) condensations and 
the hadron-quark deconfinement transition at high-density matter or the 
liquid - gas transition at subnuclear density have been studied by many au- 
thors. These phase transitions may have some implications to compact star 
phenomena, and it has been expected that recent progress in observations 
might reveal such new forms of matter. 

Such phase transitions are of the first order in most cases and the 
Maxwell construction has been applied to get the equation of state (EOS) 
in phase equilibrium; there appears a separation of spatially bulk phases 
in the mixed-phase with the equal pressure. Glendenning demonstrated a 
possibility of the structured mixed phases (SMP) in such systems by invok- 
ing the proper treatment based on the Gibbs conditions ^ , where the charge 
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density as well as the baryon-number density arc inhomogeneous. Subse- 
quently, many authors have demonstrated energetic preference of SMP and 
its existence in a wide density region, disregarding effects of inhomogeneity 
of the particle configurations and/or the electric field ^. The geometrical 
structure of SMP looks like droplets, rods or slabs as in the nuclear pasta 
phase 

The Gibbs conditions require the pressure balance and the equality of 

the chemical potentials between two phases, denoted by / and //, for phase 
equilibrium ^ . ^ For a multi-component system with more than one chem- 
ical potential, as is common in neutron-star matter, we must impose the 
equality condition for each chemical potential in order to fulfill the con- 
dition of the physico-chemical equilibrium. More definitely, we, hereafter, 
consider the charge chemical potential (/Uq) and the baryon-number chemi- 
cal potential {^b) respecting two conservation laws in neutron-star matter: 
/ig — iiq and /i^ = /i^^. On the other hand, the first condition is not 
fulfilled in the Maxwell construction, since the local charge neutrality is 
implicitly imposed, while only the global charge neutrality must be satis- 
fied. When we naively apply the Gibbs conditions instead of the Maxwell 
construction, we can sec that there appears SMP in a wide density region 
and there is no constant-pressure region in EOS. 

SMP, if exists, may have phenomenological implications on compact 
stars through e.g., glitches, neutrino opacity, gamma-ray burst or mass of 
hybrid stars. 

In this talk we address a controversial issue about the relevance of SMP, 
by taking the hadron-quark deconfinement transition as an example ^. We 
shall see that the Debye screening effects greatly modify the mechanical 
stability of SMP. In the absence of SMP we effectively recover the picture 
of phase equilibrium given by the Maxwell construction where two bulk 
phases are separated without spoiling the Gibbs conditions. 

2. Bulk calculations and flnite-size effects 

Consider SMP consisting of two phases I and II, where we assume spherical 
droplets of phase I with the radius R to be embedded in the matter of phase 
II and two phases are clearly separated by sharp boundaries. We divide the 
whole space into the eqiiivalent Wigner-Seitz cells with the radius Rw (see 
Fig.l). The volume of the cell is Vw = 47ri?^/3 and that of the droplet is 
V = 4nR^/3. 



"^We consider here matter at zero temperature 
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Figure 1. Equal droplets of the phase I embedded in the phase II (right panel), and the 
geometrical structure of the Wigner-Seitz cell (left panel) 



A bulk calculation proceeds as follows ^. For a given volume fraction 
factor / = (R/Rw)^, the total energy E may be written as the sum of the 
volume energy Ey, the Coulomb energy Ec and the surface energy Es, 

E^Ev + Ec+Es. (1) 

We further assume, for simplicity, that baryon number {p%) and charge 
(Pq) densities are uniform in each phase a, a = I, II. Then, Ey can 
be written as Ey/Vw — I^^Pb) + (1 ~ D^^Hpb) terms of the energy 
densities e", a = I, II. The surface energy £'5 may be represented as 
Es/Vw = / X 47r(T/i? in terms of the surface tension cr. The Coulomb 
energy Ec is given by 



Ec/Vw^fx^-^{pl,-pijfR'. 
The optimal value of Rd is determined by the minimum condition. 



diE/Vw) 



dR 



-0, 



(2) 



(3) 



for a given / (see Fig. 2). Since Ey does not depend on i?, we can always 
find a minimum as a result of the competition between the Coulomb and 
the surface energies, satisfying the well-known relation, Es = 2Ec. 

However, such bulk calculations have been proved to be too crude for 
the discussions of SMP. Instead, a careful consideration of the interface of 
two phases is required. As a defect of the bulk calculations they ignore the 
finite size effects. In particular, they have the inconsistent treatment of the 
Coulomb potential; they do not use the Poisson equation, so that the charge 
density profiles are assumed ab initio to be constants and the Coulomb 
potential is assumed to be 1/r. If one properly solves the Poisson equation, 
one should have the screening effect as a result of the rearrangement of the 
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Figure 2. Schematic view of the total energy and each contribution in the bulk calcu- 
lations (solid curves). Screening effect reduces the Coulomb energy, shown by the thick 
arrow. 

charge-density distribution. Hence, the radius Rd should be not too large, 
compared with the Debye screening length A^^ = X^il-^ri)"^' 



in order the above treatment to be justified, the sufhx i runs over the 
particle species. Otherwise, the Coulomb energy is reduced by the screening 
effect, which should lead to a mechanical instability of SMP in some cases 
(Fig. 2). In the case of the hadron-quark deconfinement transition, Xjj ~ 
5fm and , Afj are of the same order as A|j , for a typical density with 
IJ,B — IGeV. We shall see in the following that Rd is typically of the same 
order as Xd ~ A|), and the mechanical stability of the droplet is much 
affected by the screening effect. 

3. Mechanical instability of the geometrical structure of 



3.1. Thermodynamic potential for hadron - quark 
deconfinement transition 

In the following we consider thermodynamics for non-uniform systems. The 
situation is the same as described in Fig. 1: the phase I in the domain 
consists of M, d, s quarks and electrons and the phase II in the domain D^^ 
neutrons, protons and electrons. These phases should be clearly separated 



1/Ag - 47rQ, 



dpch 



(4) 



SMP 
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by the narrow boundary layer Ds with the width ~ ds < 1 fm due to the 
non-perturbative effect of QCD. We treat such narrow boundary as the 
sharp one (dD) with the surface tension parameter ctqcd by using the bag 
model picture, while the value of ctqcd is poorly known. We shall see that 
the Debyc screening length \]j is much longer than dg and thereby the 
introduction of the sharp boundary should be reasonable ^. 

Then, the thermodynamic potential per cell is given by a density func- 
tional ^, 

n = E[p] -p\[ dfp\ - pf [ drpf, (5) 
where E\p\ is the energy of the cell and consists of four contributions: 
E\p]= j dfei,,+,^,[pl]+ [ dfe^^„+,tJpf]+47riiVQCD+i^y. (6) 

The first two terms are given by the kinetic and strong interaction energies, 
and the Coulomb interaction energy Ey is expressed in terms of particle 
densities, 

Ey = U dfdr'^^^^^S^^, (7) 
2 J \r — r' \ 

with Qi being the particle charge (Q = — e < for the electron). 
The equations of motion are given by dSl/Spf = and written as 

= ^^1^ - '"^"(r), 7Vf'"=Qr/e, (8) 
with the electric potential V"{r): 

Thus chemical potentials pf for charged particles have values depending 
on the electric state of the phase as well as on its chemical composition. 
Actually it is sometimes called the electro- chemical potential to stress this 
fact ^ 



This treatment is also similar to the Gibbs geometrical surface ' 
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3.2. Gauge invariance 

The thermodynamic potential enjoys the invariance under a gauge trans- 
formation, V{r) V{r) — V° and tJ-i ^ fJ-f + '"V^ , with an arbitrary 
constant . Hence the chemical potential /if acquires physical meaning 
only after gauge fixing . 

Here we reconsider the Gibbs conditions and the Maxwell construction. 
As has been mentioned, on the first glance the Maxwell construction looks 
as contradicting the Gibbs conditions, especially the equilibrium condition 
for the charge chemical potential fJ-Q{= fJ-e) in our context. However, cor- 
rectly speaking, when we say /if ^ /ig^ within the Maxwell construction, it 
means nothing but the difference in the electron number density rie in two 
phases, ^ n^^; this is because rie — /^^/(Svr^), if the Coulomb potential 
is absent. Once the Coulomb potential is taken into account, using eq. (8), 
Tie can be written as 

Thus we may have /if = /if ^ and nf ^ nf ^ simultaneously, with the different 
values ofV, ^ V" (see Fig. 3). 
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Figure 3. Relation between the charge chemical potential fJ.Q{= fie) and the electron 
number density rte in the presence of the Coulomb potential V. Fulfilling the Gibbs 
conditions, fig = fJ,g, f^l = fii^ , we can change Tig in two phases as in the Maxwell 
construction, if V changes from one phase to another. 



Applying Laplacian (A) to the l.h.s. of eq. (9) we recover the Poisson 



"^Note that y = is a conventional choice in the usual treatment of uniform matter, 
while any constant is possible there. 
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equation (re D"), 

AV"{r) = 47reV'^^'"(r) = ineQfpfir). (11) 

The charge density p'^^'"(r) as a function of V"{'r) is determined by tlie 
equations of motion (8). Thus eq. (11) is a nonlinear differential equation 
for y"(r). The boundary conditions are 

yi = F", vy' = Vy", fGdD, (12) 

where we have neglected a small contribution of the surface charge ac- 

cuimnlated at the interface of the phases. We also impose the condition, 
= 0, at the boundary of the Wigner-Seitz cell, which implies that 
each cell must be charge neutral. 

3.3. Results 

The Debye screening parameter is introduced by the Poisson equation, if 
one expands the charge density in 5V" (r) = {r) — V^^ around a reference 
value V^gf. Then eq. (11) renders 

AdVif) = 47reV'^'"(V""(r) = Kef) + {^'^iV'^ir) = V^'li)f6V"{r) + 

(13) 

with the Debye screening parameter, 



dV 



(14) 

Then we calculate contribution to the thermodynamic potential of the cell 
up to 0{5V'^{r)Y . The "electric field energy" of the cell (7) can be written 
by way of the Poisson equation (13) as 

that is, in the case of imscreened approximations, usually called the 
Coulomb energy. Besides the terms given by (15), there are another contri- 
butions arising from effects associated with the inhomogeneity of the electric 
potential profile, through implicit dependence of the particle densities on 
yi'ii(r). We will call them "correlation terms", uj'^^^ — ej^jj^^^^.^ — nfpf- 

We obtain the corresponding correlation contribution to the thermody- 
namic potential Ocor = Jjji dful^^ + Jjju dfw^J^^: 

<or = ekin+str(P?(Kef)) ' ^^^^(Kef) ' P'^^'" (Kef )^ref 

V;gfAV-"(^") , i'^''iK"i)?iSV"{r^? , .... 
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where we also used eqs. (13) and (14). In general V^^f ^ V^^f and they may 
depend on the droplet size. Their proper choice should provide appropriate 



convergence of the above expansion in 5V{r). Taking V^^f 
const we find 



except an irrelevant constant. 



(17) 




Figure 4. Dimensionless value of the thermo- 
dynamic potential per droplet volume. Solid 
lines are given for / = 0.5 and dashed lines for 
/ = 1/100. The ratio of the screening lengths of 
two phases, ao = Aj^/Aj^, is fixed as one. ^ is a 
dimensionless radius of the droplet, § = R/Xj^, 
with Aj-, ~ 5 fm in this calculation. See text for 
further details. 



For given baryon-number 
chemical potential /is and 
charge chemical potential /ig, 
all the particle chemical po- 
tentials fii can be represented 
in terms of /is and /xq(= /ie) 
with the help of the chemical 
equilibrium conditions: 

- fJ-s + yUe = 0, Ms) 
^^n ^ ^J.p + tJ-e, (18) 

in each phase and 

= ^„ = 2^id + ^J'u, (19) 

at the boundary. 
Then particle number den- 
sities Pi are represented as 
functions of /zb, /zq and the 
Coulomb potential V, due to 
the equations of motion. Sub- 
stituting them in the Poisson 
equation (11), we can solve 
it with the proper bound- 
ary conditions (12); note that 
Vy = at the boundary of 
the Wigner-Seitz cell provides 
us with another relation be- 
tween II B and pq. 



Thus we eventually have the density profiles of all the particles for 
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given density or the baryon- number chemical potential Hb- In Fig- 4 we 
demonstrate the radius (R) dependence of the total thermodynamic po- 
tential per droplet volume for the case of spherical droplets, JStot/Z^o = 
(cy + + wJoj. + wJoj + ?s) /Po, given by the sum of partial contributions, 
where tilde denotes each quantity scaled by the droplet volume V = 47ri?^/3 
and is a typical quantity with the dimension of the energy density^. 
Preparing some wide range for the value of the surface tension parameter 
o'QCD, Pi oc (TQCD^, we present two cases of /, / = 0.01, 0.5. 

The label "C" is given for reference to show the previous non- 
selfconsistent case, where the Coulomb potential is not screened see Fig. 2. 
We can see that only in the limit of / <C 1 and <C A^, we are able 
to recover this case. The "c.m." curve shows the partial contribution to 
the thermodynamic potential, Ce.m./Po = (^v +^s)/Poi ignoring correlation 
terms. Comparing these curves we can see how the screening effect changes 
the thermodynamic potential: wc can see that the minima at the "c.m." 
curves disappear already at f3i > 0.03, corresponding to unphysically small 
o'QCD ~ several MeV-fm~^. However, the correlation energy gives a sizable 
contribution to allow the minimum for larger value of (Jqcd- Consequently, 
the minimum totally disappears between /3i = 0.1 and /?i = 0.5, which may 
be interpreted as 10 < cfqcd < 50(MeV-fm~^) in this calculation. Thus 
we have seen a mechanical instability of the droplet for the medium values 
of (JQCD, which might be in the physically meaningfuU range. 

4. Summary and Concluding remarks 

In this talk we addressed an issue about SMP at the first order phase 

transitions in multicomponent systems, like in neutron-star matter. We 
have studied a so called "contradiction" between the Gibbs conditions and 
the Maxwell construction extensively discussed in previous works. We have 
demonstrated that this contradiction is rc;solvcd if one correctly takes into 
account the difference in the "meaning" of the chemical potentials used in 
the two approaches: the different values of the electron chemical potentials 
in the Maxwell construction and the ones used in the Gibbs conditions do 
not contradict each other if one properly takes into account the electric 
field. 

We have presented a framework based on the density functional theory 
to describe thermodynamics in the non-uniform systems. The Coulomb 
potential is properly included and particle density profiles are consistently 
determined with the Poisson equation. 
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Taking the hadron-quark deconfincmont transition in high-density mat- 
ter as an example, we have demonstrated the importance of the Debye 
screening effect, which is a consequence of the above treatment. With a 
numerical example, we have seen that the screening effect gives rise to a 
mechanical instability for realistic values of the surface tension parameter 
of (JQCD- In this case we may effectively recover the picture given by the 
Maxwell construction, where the phase separation of two bulk phases arises. 

Our framework is rather general and it may be applicable to any first 
order phase transition, e.g. the liquid-gas phase transition at subnuclear 
density ^. 
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